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Abstract 

Motivated for the fault tolerant quantum computation, a quantum gate by adi- 
abatic geometric phase shift is extensively investigated recently. Due to the 
enviromental noise, the decoherence time for a quantum state can be very short. 
In the adiabatic quantum computation, if we change the Hamiltonian slowly, 
the decoherence error is increased, if we change the Hamiltonian faster, the 
adiabatic condition is broken and the error caused by the nonadiabaticity is in- 
creased. In this paper, we give a nonadiabatic scheme for the geometric phase 
shift and conditional geometric phase shift with the nuclear magnetic resonance 
system (NMR). Essentially, the new scheme is simply to add an appropriate ad- 
ditional field. With this additional field, the state evolution can be controlled ex- 
actly on a dynamical phase free path. Geometric quantum gates for single qubit 
and the controlled NOT gate for two qubits are given. In our non-adiabatic 
scheme, the geometric quantum gates can run as fast as the usual quantum 
gates. 
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I. INTRODUCTION 



It has been shown that, a quantum computer, if available, can perform certain tasks much 
more efficiently than a classical Turing machine. A fault-tolerant quantum logic gate [|IJ is the 
central issue in realizing the basic constituents of a quantum information processor. Quantum 
computation via the controlled geometric phase shift 0] provides a nice scenario to this purpose. 
Due to its geometric property, a geometric phase |||| shift can be robust with respect to certain 
types of operational errors. In particular, suppose that the spin(qubit) undergoes a random 
fluctuation about its path in the evolution. The final value of the geometric phase shift will 
not be affected provided that the random fluctuation does not change the total area on Bloch 
sphere @. 

Recently, it was reported @J^-0 that the conditional Berry phase( adiabatic cyclic geometric 
phase) shift can be used in quantum computation. In particular, in ref |J, an experiment was 



done with nuclear magnetic resonance(NMR) JlO-14] under the adiabatic condition. 

Demand on both running speed and precision of every gate in a quantum computer is quite 
high. By the available technology at present ( and at near future), due to the enviromental 
noise, distortion to a quantum state increases fast with time. This requires all quantum gates 
complete the task as fast as possible. However, the schemes raised previously required the 
adiabatic condition. That is to say, if we run the gate too fast (i.e. change the Hamiltonian 
rapidly), the distortion from the nonadiabaticity must be significant. It has been reported in 
the recent experiment that 0, due to the requirement of adiabatic condition, faster running 
speed causes severe distortions to the results . The adiabatic condition makes the fast speed 
and high precision conflict each other in the conditional Berry phase shift gate. Increasing the 
running speed of the logic gate can exponentially increase the power of a quantum computation 
such as quantum factorization fT5| . It seems that the geometric phase shift gate will only be 
practical if it can run at a speed comparable to that of the usual quantum gate. Therefore 
one is tempted to find an easy way to make the conditional geometrical phase shift through 
nonadiabatic state evolution. In this paper, we present an easy scheme to make the geometric 
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phase shift nonadiabatically with NMR. In section II we study the state evolution on a cone. 
We give a simple way to nonadiabatically control the state evolution on the cone. With an 
appropriate setting, the state evoleves on a dynamical phase free path. In section III we im- 
plement the result of section II to the NMR system. A scheme for nonadiabatic conditional 
geometric shift is given. In section IV we investigate the nonadiabatic geometric realization of 
various types of quantum logic gates, including the phase shift gate, the Hadamad gate and 
the NOT gate for single qubit, and the C-NOT( controlled-NOT) gate for two qubits. We give 
a concluding remark in section V. 



II. EXACT STATE EVOLUTION ON THE CONE. 

It is well known that a spin half nucleus can gain a geometric phase shift in the conical 
evolution. We now assume that the initial spin state is on the cone and we will demonstrate a 
method to make it evolve nonadiabatically. The Hamiltonian for a spin in a constantly rotating 
magnetic field Q is 



H(t) = [u a z + ui<r x (t)]/2. 



(1) 



Here u x (t) 



e-^ ^ 



V 



ljo is the amplitude of vertical field and U\ is the amplitude of the 



horizontal field, which is rotating around z— axis in the constant angular speed 7. The initial 



state iV'o) is an eigenstate of Hq = H{0)- Explicitly, 



\ip ) = cos^| t) + sin 1 1 |) 



(2) 



and 



cos 6 = u / \J(jJq+ oof. 



(3) 



1 We shall extend it to the case of arbitrary time-dependent rotating speed in the end of this section. 
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To know the cause of the state distortion, we study the time evolution first. (The similar idea has 
been used for the adiabatic rotational splitting with NQR previously (JTjJ.) The time evolution 
operator U generated by the Hamiltonian H is determined by the time-dependent Schrodinger 
equation 

i^U = H(t)U. (4) 

Solving this equation we obtain 

U(t) = e -^t e -i(H"-~f^/2)\ (5) 

In particular, at time t = t = 27i~/7, when the external field completes a 2ir rotation, the state 
is evolved to 

^(r) = e - m - i{H °-^ /2)T \i) ). (6) 

The adiabatic approximation assumes that at time r, the state completes a cyclic evolution 
provided that 7 is small. However, due to the decoherence time limitation, 7 cannot be too 
small. The non-zero 7 distorts the state in the evolution and makes it noncyclic at time r . This 
non-zero 7 in the adiabatic approximation causes two types of errors. One is that the noncyclic 
completion of the evolution at time r will cause further errors in the succeeding operations in 
the sequence. (A sequence of operations is used to remove the dynamic phase in realizing the 
geometric quantum gate ||.) The other is that the geometric phase acquired over period r is 
not — 7r(l — cos#), as for the ideal adiabatic cyclic evolution. 

Now we give an easy way to exactly control the state evolution on the cone. Here the 
external field can be rotated around z— axis arbitrarily fast. We use l^o), the eigenstate of H 
for the initial state. We switch on a static vertical magnetic field u z = 7 while the external 
field is rotated around z— axis. In this way, the new time-dependent Hamiltonian is 

H w (t) = ^[(uj + j)a z + cu 1( T x (t)} . (7) 

The time evolution operator Uw(t) generated by this Hamiltonian satisfies the Schrodinger 
equation 



i^-U w (t) = H w (t)U w (t) (8) 

with the boundary condition U w (0) = 1. The state at time t is related to the state at time 
by \ip(t)) = Uw(t)\i>o), 4>o = ^(0). Denoting R = e n ' CTz//2 we obtain the following equivalent 
equation 

id(RU w )/dt = H (RU W ). (9) 
Now that H is time-independent we have 

RU w {t) = e~ iHat . (10) 

This is equivalent to 

U w (t) = e -*7t*»/2 e -itfot, (11) 
Consequently at any time t, state 

\^{t)) = e- a V^/Vo) (12) 

is exactly the instantaneous eigenstate of Hamiltonian H(t), where A is the eigenvalue of H 
for eigenstate |^ ), 

A = ±yju$ + uj\. (13) 

In particular, at time r, 

|^(r)) = e— ^|^ ). (14) 

It only differs to |^o) by a phase factor. Note that the time evolution operator Uw(t) here is 
generated by the Hamiltonian H w (t) instead of H(t). 

Here the additional vertical field 7 plays an important role. Without this field, the time 
evolution operator generated by H{t) is 

U(t) = e -*t*./2 e -iHit ( 15 ) 



and Hi = H — ja z /2. Obviously, this time evolution operator U(t) will distort a qubit with 
the initial state \ipo) in the evolution. But this U(t) can exactly control the qubit with the 
initial state the eigenstate of Hi. If initially we set the spin state to \ipi), then the spin 
will evolve exactly on its cone without the additional field 7. So, we have two ways to control 
the spin evolution exactly. We can use the additional field 7, if the initial state is set to be 
|^o) ■ Alternatively, we can also set the initial state to be \ipi) and then we need not add any 
additional field when the external field is rotated. In this letter, we adopt the former one, i.e. 
we set the initial spin state to \ipo). 

The additional magnetic field 7 can make the state \ipo > keep up the rotating field exactly in 
the evolution. So the state completes a cyclic evolution when the field completes a 2tt rotation. 
However, after a cyclic evolution, the total phase shift for the state in general includes both 
dynamic and geometric contribution |I7|| . It is possible to restrict the state evolution on the 
dynamic phase path if we choose certain specific value of the rotating speed 7( and also the 
additional field 7). For this purpose, we require 



7 cos 6 = — yco>o + w i- (16) 

The solution is 



7=-^°. (17) 
^0 



The negative sign in the right hand side indicates that the additional field 7 is anti-parallel 
with the field uq. By this setting, the total magnetic field in Hyy is always "perpendicular" to 
the state vector expressed in the Bloch ball. Details of this are shown in Fig. 1. One can easily 
show that the instantaneous dynamical phase 

< MUl{t)H w {t)U w {t)\^ >= 0. (18) 

We have assumed above that the field rotates at a constant speed. Actually, we can easily 
modify the above scheme for arbitrary time-dependent rotating speed 7 a (£), Jo 1a{t)dt = 2n. 
In this case, we need only change the term jt in Hw(t), Uw(t), R and ip(t) into J t 'j a ( y t')dt' 



accordingly. Consequently, the the additional vertical field is now a time- dependent field uj z (t) = 
7 a (t) instead of a static field. To remove the dynamical phase, we can use the time-dependent 
magnetic field coi(t) and u> (t). We require 

LOx[l) 

and 

wi(t)/w (*) = wi(0)/w (0) = tanfl (20) 

This extension to the time- dependent case can be important in case it is difficult to rotate 
the field (or the fictitious field ||) in a constant angular speed. Punctual results can be obtained 
here through the exact feedback system where the value of additional vertical field is always 
instantaneously equal to the angular velocity of the rotating field. 

Thus we see, by adding an additional magnetic field that is equal to the rotating frequency 
of the external field, we do get the exact result. Obviously if we rotate the field inversely, 
additional vertical field should be in the inverse direction(— z) accordingly. 

III. THE NONADIABATIC CONDITIONAL GEOMETRIC PHASE SHIFT. 
A. NMR system and the rotational framework. 

Consider the interacting nucleus spin pair (spin a and spin b) in the NMR quantum compu- 



tation |iOHT3|. If there is no horizontal field the Hamiltonian for the two qubit system is 



Hi = -{u a a za + u b <j zb + Ja za ■ a zb ) (21) 
, where u a (uj b ) is the resonance frequency for spin a(b) in a very strong static magnetic field(e.g. 



u a can be 500MHz [Q]) , J is the interacting constant between nuclei and a za = a zb = a z = 
. After adding a circularly polarized RF field in horizontal plane, the Hamiltonian 
for spin a in static framework is 



d o \ 



V " 1 / 



1 



1 , 



1 



H' = -u cr z + -J a a z + -ui 



e 



-iu)'t 



\ 



(22) 



/ 



Here u' a and W\ are the angular frequency and amplitude of the RF field respectively, 



= u a - J a ± J 



(23) 



and the " ± " sign in front of J is dependent on the specific state of spin b, up or down 
respectively. To selectively manipulate spin a, we shall use the rotational framework that is 
rotating around z— axis in speed u' a . We assume ui' a close to u a but quite different from uib- The 
Hamiltonian for spin a in the rotational framework is 



1 



1 



H a = R'H'R + i{dR' /dt)R'- 1 = -u a z + -u l( r x {t) = H{t) 



(24) 



and R! = e lu) ' a(Jzt / 2 . Note here H is dependent on state of spin b through ujq. In the rotational 
framework, if the horizontal field is rotated in the angular speed 7, the Hamiltonian for spin a 
is just H(t), as defined in eq(|l]). In the NMR system, we require j^o) to be the eigenstate of H 
in rotational the framework. Previously 0, state l^o) f° r spi n a was produced adiabatically. 
In the following we demonstrate how to nonadiabatically produce the state l^o). 



B. Creating the conditional initial state with NMR. 

In rotating the state around z— axis, we have started from state \ipo), which is the eigenstate 
of H . State j^o) can be created from the spin state | |) or | j). We have a nonadiabatic way to 
create state |f/>o) for spin a. Note that we are working in the rotational framework. We denote 
5 = uj a — u' a . We can use the following sequence of operations to create the conditional angle 
#.(For simplicity we will call the following sequence as S operation later on.) We have 



S 



(25) 



Here all terms inside the [■ ■ •] represent the Bloch sphere rotation angles caused by RF pulses. 
The superscripts indicate the axis the Bloch sphere is rotated around. J'((p±(t c )) is the time 
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evolution over period t c by the Hamiltonian |(5 ± J)a z . This evolution rotates spin a around 
z— axis for an angle ip± — (5 ± J)t c . (Or </? ± + 7r, if spin a is down initially. For clarity we omit 
this case and alway assume spin a is initially up.) After this S operation, the angle between 
spin a and the z-axis is 9± = f — (<// + <^±)(see Fig. 2). To ensure the state to be eigenstate of 
Hamiltonian H after S operation we require tan((/?' + Jt c ) = and tan(<// — Jt c /2) = 
simultaneously. This is equivalent to 

Jt c = (arctan ^ - arctan ^) 12 



(26) 

= (arctan + arctan ^) /2 
From this constraint, given specific values of 5, J and u>i, we can easily obtain the scaled time 
J ■ t c (Fig. 3) and the angle (p'(Fig. 4) in controlling the S operation. In particular, for 
5/J= 1.058, the value adopted in the recent NMR experiment, the relation curves for Jt c vs 
uj 1 and ip' vs uj\ are shown in Fig. 5. 

C. The nonadiabatic conditional geometric phase shift. 

After the S operation, the state of spin a is an eigenstate of Hamiltonian Hq no matter 
whether spin b is up or down. After this S operation we can use the additional magnetic 
field to nonadiabatically control the state evolution on the cone. However, here we need again 
remove the dynamical phase contribution. Note the situation here is different from the single 
qubit case. Qubit b could be either up or down. We need choose the appropriate 7 value so 
that the instantaneous dynamic phase for qubit a is always zero no matter qubit b is up or 
down. So in this case the following simultaneous conditions are required. 

7 cos^± = —\J [5 i J) 2 + L)j\. (27) 

Setting 



Ul = V5 2 - J 2 . (28) 

will make the simultaneous conditions hold. With this setting we propose the following scheme 
to make the nonadiabatic conditional geometric phase shift. 
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s 



H 



(29) 



The term 



( 7 \ 



represents doing operation C with an additional vertical field 7. C represents 



rotating the external field around z— axis for 2n in a uniform speed 7. We should choose the 
rotating direction so that additional field 7 is anti-parallel the field 5± J, (5 > J). Since qubit 
a evolves on the dynamical phase free path, the scheme raised here can remove the dynamic 
phase and retain only the geometric phase. It can be shown that the geometric phase acquired 
after the operation is T + , — r + , T_ and — T_ respectively for the four different initial state 

(in>,lu>,UT>,lu», 



r± = — 7T — 7T COS#± 



(30) 



and 



COS 04 



5±J 



I5±J 
25 



(31) 



IV. GEOMETRIC QUANTUM LOGIC GATES 

A. Geometric gates for single qubit 

To make a quantum computation, we need both single qubit gates and the two qubit C- 
NOT gate. It is impossible to make any single qubit geometric gate with NMR in the adiabatic 
condition because the interaction among different qubits cannot be switched off. Though the 
interaction is small, the effects from the interaction may accumulate significantly after a long 
term. Strictly speaking, the adiabatic geometric quantum computation raised previously @ is 
only a partial geometric quantum computation. However, in the our non-adiabatic scheme, we 
can easily make single qubit gate. In this nonadiabatic case, the effect of the interaction among 
qubits is negnigible because the operation time is very short. 
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In section II, we have required the initial state of the qubit be the eigenstate of the initial 
Hamiltonian, H Q . However, this requirement is not necessary in making the geometric quantum 
gates, instead of making geometric phase shift to the qubit directly. If we are only interested 
to the geometric quantum gates instead of making the phase shift to the qubit itself, we can 
directly use arbitrary initial state. We now demonstrate this by some examples. Without any 
loss of generality, we suppose the initial state of the single qubit is | |> or | |>. We start from 
the phase shift gate. We set the initial magnetic field uj in the z direction. We then rotate uj 
anti-clockwise around z' axis. The rotating angular velocity is 7. The angle between z' and 
z axis is #o( see Fig. 6). Meanwhile we add the additional field 7. The rotating speed and 
additional field should satisfy eq. (|T6|) . Note that in this case, the additional field is in the —z' 
direction. After a 2tt angle rotation the initial states will gain phase shifts according to the 
following equation 



(32) 







= W p (6 ) 






V 


l>) 




\\i>) 



and 



P V u o) 



-2m cos t 



\ 









3 27ri cos 9q 



(33) 



The phase shift here is determined by area enclosed by the state evolution loop on Bloch 
sphere, i.e. by angle 6q. Since #0 can take arbitrary value, we can make arbitrary phase shift 
gate through the above equation. 

Now we consider the Hadamad gate, which is another fundamentally important gate in 
quantum computation. We come back to z axis. Now 8q is the angle between magnetic field uj 
and z axis. With the constraint of eq([TJ]), suppose the additional magnetic field and rotating 
speed of uj is 7. State | T> or | [> can be written in the following linear superposed forms 
respectively 



* $0 1 . $0 1 

T>= cos — 1+ > — sin — |— > 



(34) 
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and 



l>— sm—\+ > +cos — |- > 



(35) 



where the states |± > are the spin states parallel or anti-parallel to the direction of u = 
\JujI + uo\. After the field rotates one loop around z axis, the initial state | |> or | |> changes 
by the following formula 



V U>, 



V ll>y 



(36) 



where 



cos y -sin-^ 1 



sin cos 



2 J 



e 



cos y sin y 



sinf cosf j 



(37) 



i.e. 



W s = 



( 

V 



/ 



(38) 



cos r cos 9 + i sin T sin # i sin T sin # 

i sin T sin # c °s r — i sin 7 sin # 

Here T is the geometric phase gained by state |+ > after it completes the loop around z axis. 
If we take | sinTsin^ol = ^(root ^0 for this equation obviously exists), this W s is 



W. 



V2 



sO 



e i<t> I 



\ 



(39) 



Here = tan 1 smrsi ° 6>0 . Combining this W s0 gate with an appropriate phsae shift gate we can 
obtain the Hadamad gate. Explicitly we have 



/ e -»(*- w /2) 0^ 



W. 





_ V2 


(1 


1 \ 


v l ) 


2 




-v 



= w H . 



(40) 



Obviously, by combining the gate Wh and the W p , we can also make the NOT gate. It is 



w„w r (e„ = cos- 1 i)w„ = . 



(41) 



V 1 °/ 
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The NOT gate of single qubit played an important role in the adiabatic scheme to cancel the 
dynamical phase 0|§. 

So far we have demonstrated the geometric realization for several types of logic gates in 
single qubit operation. Here geometric phase is the only contribution in the time evolution. 
Actually, all single qubit gate can be realized in a similar way. The running speed of the above 
geometric gates can be in principle arbitrarily fast. 



B. Geometric C-NOT gate 

With the conditional geometric phase shift, it is very easy to construct the C-NOT gate. 



To do so, we take S = ^-J, then we have 



7T 

2 



(42) 



and 



S 





( 




s- l = 




\ 



e i(r_-7r/2) o 
Q e -i(r_-7r/2) 



(43) 



if qubit b is up. And 



S 



s 



-1 



( e ir 



V 



o 

e~ ir - 



\ 
) 



(44) 



if qubit b is down. We can use a specific phase shift gate W p on qubit a to convert the above 
form of conditional phase shift gate into 

/ 



S 



/ 7 \ 



s- 1 



e~ ir - 














e ir - 














e~ ir - 

















i 00 
10 
1 



(45) 



Sandwitching it by two Hadamad gates to the single qubit a we obtain the C-NOT gate, i.e. 



13 



V2 
2 



^110 0^ 



1-10 
11 
1-1 



i 
10 
1 



x/2 
2 



^1 10 0^ 



1-10 
11 
1-1 



^ -i 0^ 



-i 
10 
1 



(46) 



V. CONCLUDING REMARK. 

To make the conditional geometric phase shift gate we need be able to control two elementary 
operations. One is to exactly control the cyclic state evolution on a cone; the other is to produce 
the initial state on the cone with an angle 8± conditional on the other bit. Previously [[|, these 
two elementary operations were done adiabatically. We have shown that both of these two tasks 
can be done nonadiabatically. We have demonstrated the non-adiabatic geometric phase shift 
and quantum logic gates for both single qubit and two qubit case. The single qubit geometric 
gate can run arbitrarily fast in princilpe. The running speed of two qubit gate is is only limited 
by the J, the value of the qubits interraction. This is the case for all normal conditional quantum 
gates. So it possible to run the geometric quattum gate in a speed comparable to that of the 
normal quantum gate. The idea on nonadiabatic geometric phase shift gate demonstrated by 
the NMR system here should in principle also work for the other two level systems, such as 
Josephson Junction system || and the harmonic oscillator system [§,0], where the decoherence 



time can be much shorter. It should be interesting to study the nonadiabatic scheme for the 
nonabelian geometric quantum computation [[H5|-|2"2" 
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FIGURES 




yV/ 

FIG. 1. Control the conical evolution exactly. The initial state is the eigen state of Hq. 
uj is the field in Hq. uj is rotating around z axis in the direction indicated by the arrow. Additional 
field 7 is equal to the rotating speed of uj. With this additional field, the state will evolve keep up the 
field rotating uj exactly. We can see that, choosing an appropriate 7, the total field is instattaneously 
perpendicular to the field u>, thus eigen state of Hq evolves in a dynamic phase free path. 
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[7t/2] y — /'(qKtc))— [-8tc]^ 




a b c 



[Jt/2]U- [-(p'] y -^ 



G /I 




cf e f 

FIG. 2. Using 5 operation to produce the eigenstate of Hq nonadiabatically. Picture 

a shows the possible initial state for spin a. The up(down) arrow on the solid line represents the 
up(down) state for spin a while spin b is up. The up(down) arrow on the dashed line represents the 
up(down) state for spin a while spin b is down. In picture f, angle GO A is 6+ = arctan j^j and angle 
GOB is 6L = arctan ^j. 




FIG. 3. The time control in S operation. Parameter ui/J varies from to 10, 5/ J varies 
from to 5. Vertical axis is for the scaled time Jt c . 
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FIG. 4. The rotating angle control in S operation. Parameter uj±/J varies from to 10, 5/ J 
varies from to 5. Vertical axis is for the angle (p'. 



1.4 - 

1.2 \ 

1 4 

0.8 - ' 




1 1 1 1 1 1 1 

1 2 3 4 5 

FIG. 5. The time and rotating angle control in S operation for specific 6/ J value. 

Horizontal axis represents uj\/J. Vertical axis is for the scaled control Jt c (the solid line) or the angle 
(//(dashed line). Here 5/ J = 1.058, as for the experimental condition[2]. By this figure, given the 
specific values of lo±/J, we can always find the corresponding point in the two curves thus we can take 
the suitable time control(for t c ) and rotation control(for ip f ) in the S operation. 
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FIG. 6. Rotating the field around z' axis. Here the field u> is in +z direction, rotates around 
z 1 axis. Additional field 7 is in —z 1 direction. 
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